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O ■ Abstract 

The hadronic particle production data from relativistic nuclear Pb-Pb 158 A GeV collisions are 
f""**. ' successfully described within the chemical non-equilibrium model, provided that the analysis treats 

■ SI and Q abundances with care. We further show that there is a subtle influence of the Coulomb 

potential on strange quarks in quark matter which is also seen in our data analysis, and this Coulomb 
effect confirms the finding made by chemical analysis in the S-Au/W/Pb 200 A GeV collisions that 
' the hadron particle source is deconfined with respect to strange quark propagation. Physical freeze- 

t^J- ■ out conditions (pressure, specific energy, entropy, and strangeness) are evaluated and considerable 

C*~) ' universality of hadron freeze-out between the two different collision systems is established. 
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' 1 Introduction 

Oh! 

Qh! Intense experimental and theoretical work proceeds to explore the mechanisms of quark confinement 
effect and the properties of the vacuum state of quantum-chromodynamics (QCD), the non-Abelian 
gauge theory of 'color' charges Q. Relativistic energy nuclear collisions are the novel experimental tool 
developed in the past decade to form, study, and explore the 'melted' space-time domain, where we hope 
to find, beyond the Hagedorn temperature Th c± 160 MeV ||, freely propagating quarks and gluons 
in the (color charge) plasma (QGP). Analysis of hadronic particle production in relativistic nuclear 
collisions offers an opportunity to explore the mechanisms of quark confinement at the time of final 
state QGP freeze-out. We present here a progress report of such ongoing analysis of the recent and 
forthcoming experimental study of the Pb-Pb collisions at 158 A GeV Q, which follows on our now 
complete analysis of the S-Au/W/Pb collision systems Q. 

There is little doubt that in the early Universe QGP was the transient state of matter, and that 
only about 20-40 //sec into evolution did our present confining vacuum freeze-out from the primordial 
QGP-form. The issue is if in laboratory experiments we can indeed form and study the primordial QGP 
phase H|. In some aspects, such as specific entropy and baryon content, notable differences between the 
laboratory QGP state and the early Universe conditions are present. Also, the small, nuclear size of the 
nuclear collision 'micro'-bang implies a short laboratory lifespan r q ~ 0.5 • 10 _22 sec. There is also the 
difficult problem of proving the fundamental paradigm beyond a shade of doubt: can there indeed exist 
a locally deconfined space-time domain with energy density exceeding by an order of magnitude that of 
nuclear matter? Much of the current effort is solely addressing this question. Among several proposed 
approaches to search for and study the deconfinement, our work relies on the idea of strangeness flavor 
enhancement, and the associated enhancement of (strange) antibaryon formation ||. This signature 
can be combined, in the present analysis, with the study of global particle abundance which represents 
the entropy contents of the deconfined phase [0. || . 



1 



2 



J. Letessier and J. Rafelski 



We address here 15 presently available particle yield ratios obtained in central Pb-Pb 158 A GeV 
collision experiments carried out at CERN-SPS. Our analysis addresses results of experiments NA49 
|)|, [Tcfl , and WA97 p| , we have not used results from NA44 [12], being uncertain about the impact 
of the cascading weak decay contamination of the hadronic ratios, which are quite significant. Four 
WA97 data points involve O and particles and even a cursory study of these abundances suggests that 
these entirely strange particles are not falling into the same systematic class, a fact also visible in their 
unusual spectral slopes. We believe in view of the difference in systematics that it is appropriate also 
to consider the data excluding the and f2 yields from analysis. In that case our analysis contains 11 
relative experimental particle yields. However, 4 of these ratios originating in the same experiment are 
related by a simple algebraic constraint (e.g., A/ A = A/S • H/E • H/A), leaving us with ten independent 
measurements. As we shall see, there are up to 5 parameters in our description. In the different 
analysis discussed here we thus have no less than 5 independent degrees of freedom: ndof > 5. Since 
we address ratios of strange particles as well as ratios of total abundances of positive and negative 
hadrons, we combine in the present analysis strangeness observables with the entropy enhancement Q. 
Underlying our data analysis is the assumption of local thermal (i.e., energy equipartition) equilibrium. 
Both, the thermal appearance of produced particle spectra [^, [|, [l3|], and the qualitative and systematic 
agreement over many orders of magnitude between properties of a thermal hadron system and the 
experimental hadron abundance yields [14|, provide a solid foundation for the assumption of (near) 
thermal equilibrium in the dense hadronic matter. 

In our most recent theoretical approach there is a key refinement not present in earlier work Q : we 
do not assume chemical equilibrium even for light quarks. Consideration of non-equilibrium chemical 
abundance for strange quarks allowed to analyze accurately the experimental particle abundance data 
and to characterize precisely the properties of the presumably deconfined source []|, [?], 1C, |l^, 15, 16, 



17, [L8|, 19, pQ ]. The mechanisms of chemical equilibration requiring reactions which change particle 
abundances are today much better understood theoretically than those responsible for what is believed 
to be much faster thermal (kinetic) equilibration, where momentum exchange between existent particles 
is the key mechanism. It is hard to understand, why we, along with others, have maintained in the 
past in our analysis the point of view that only strangeness has the opportunity to be off-equilibrium 
in chemical abundance. Indeed, if QGP is the particle source the need to assimilate by fragmentation 
the gluon content must generate excess light quark abundance. Since we did not allow for this freedom 
in earlier data analysis, the data were (equally badly) also described by a high temperature source 
model ||]- This is not the case once also for light quarks the chemical nonequilibrium is introduced, 
only the low temperature freeze-out alternative has a convincing statistical significance 

Hadronic particles we observe are either emitted directly or are descendants of other hadronic pri- 
maries produced near or at surface of the dense matter fireball. Local rest-frame temperature T and 
local collective flow velocity v c characterize the momentum space distribution of particles emerging from 
the surface region of the fireball. In addition, each surface volume element is characterized by chemical 
abundance factors we shall discuss in more detail below. Only particles of similar mass and cross section 
experience similar drag forces arising from local flow of matter and hence ratio of their abundances in 
some limited region of phase space for not too small momenta is expected to remain unaltered by v c . 
Since the surface vector flow v c is a priori largely unknown, in order to use limited 'windows' of particle 
momenta (rapidity y and transverse mass m±) for determination of chemical freeze-out properties of the 
source, only ratios of 'mutually compatible' particles can be considered, aside of ratios of total particle 
abundances. Our recent study of the flow effect shows that this procedure is adequate, though with the 
ongoing improvement of the data sample we will soon have to include this additional freeze-out flow 
velocity parameter explicitly in the analysis pl| . 

In the following section, we shall briefly summarize the theoretical foundations of the current analysis 
and discuss the impact of Coulomb effect in QGP on the expected values of statistical parameters, and 
we shall describe the parameters that we obtain. In section ||, we will discuss the particle abundances 
and other related results, such as the physical properties of the freeze-out dense matter. We conclude 
with a few brief remarks about the relevance of our analysis for the search for QGP. 
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2 Statistical model and Coulomb effect 



dN i = ^TZW AiVidt ■ W 



The thermal production yield dNi of particles emitted within the time dt from a locally at rest surface 
element dS is: 

dSd?p 

Here Vi = dz/dt is the particle velocity normal to the surface element dS = dxdy. For a thermal quark- 
gluon gas source and allowing for recombination-fragmentation of constituents and detailed balance, the 
complete phase space occupancy factor A$ is given by: 

Ai = giKne- E ^ T , X i = l[X j , H = l[ lj, Ei = Y, Ej, (2) 

jei j& jei 

where gi is the degeneracy of the produced particle, and E{ its energy. The valance quark content {j} in 
hadron {i} is implied in Eq. The fugacities Xj arise from conservation laws, in our context, of quark 
(baryon) number and strangeness in the particle source. X q = e^ q / T is thus the fugacity of the valance 
light quarks. For a nucleon Xn = A^, and hence the baryochemical potential is: /x& = 3u, q . Similarly, for 
strange quarks we have A s = e^ a l T . For an antiparticle fugacity Xj = X^ 1 . Some papers refer in this 
context to hyper-charge fugacity As = X q /X s , thus fis = — Ms- This is a highly inconvenient historical 
definition arising from considerations of a hypothetical hadron gas phase. It hides from view important 
symmetries, such as X s — > 1 for a state in which the phase space size for strange and anti-strange quarks 
is the same: at finite baryon density the number of hyperons is always greater than the number of anti- 
hyperons and thus the requirement (N s — N s ) = can only be satisfied for some nontrivial X s (X q ) ^ 1. 
Thus even a small deviation from X s — > 1 limit must be fully understood in order to argue that the 
source is deconfined. Conversely, observation of X s ~ 1 consistently at different experimental conditions 
is a strong and convincing argument that at least the strange quarks are unbound, i.e., deconfined. 

As this discussion also illustrates, the parameter X s does not regulate the total number of s-s quark- 
pairs present in the system. More generally, any compound object comprising a particle-antiparticle pair 
is not controlled in abundance by a fugacity, since the formation of such particles does not impact the 
conservation laws. In consequence, the abundance of, e.g., neutral pions comprises no quark fugacity. 
This is the reason to introduce an additional chemical phase space occupancy factor 7j: the effective 
fugacity of quarks is X^i and antiquarks A~ 1 7j. This parameter allows to control pair abundance 
independently of other properties of the system, and in particular temperature. For 7$ — ► 1 one reaches 
a entropy maximum ||, corresponding to the 'absolute' chemical equilibrium J22|J . Therefore the factor 
7i is called the (chemical) phase space occupancy factor. 

A time dependent build up of chemical abundance was first considered in the context of microscopic 



strangeness production in QGP |2^, [23|], after it was realized that strange flavor production occurs at 
the same time scale as the collision process. More generally, one must expect, considering the time 
scales, that all quark flavors will not be able to exactly follow the rapid evolution in time of dense 
hadronic matter. Moreover, fragmentation of gluons in hadronizing QGP can contribute additional 
quark pair abundance, conveniently described by the factor 7$. It is thus to be expected that also for 
light quarks the chemical phase space occupancy factor 7 3 ^ 1. Introduction of the factor 7 g leads to a 
precise chemical description of the S-Au/W/Pb 200 A GeV collisions H, which was not possible before. 
The tacit choice j q = 1 has not allowed previously to distinguish the different reaction scenarios in 
Pb-Pb collisions ||, where we found analyzing the experimental data that hadronic particles could be 
born either at high temperature T ~ 300 MeV or at expected hadronization temperature T ~ 150 MeV. 
Introduction of 7 9 , along with improvement in precision, allowance for quantum (Bose/Fermi) corrections 
to the Boltzmann distribution functions, and a greater data sample, allowed us moreover to recognize 
the systematic difference between data points containing, and resp., not containing Q, Q, allowing us to 
develop the precise analysis here presented. 

In another refinement both u, (i-flavor fugacities X u and A^ can be introduced, allowing for up-down- 
quark asymmetry @. We recall that by definition 2/i g = ^ + fi u , thus X q = \fX u Xd. For the highly 
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Coulomb-charged fireballs formed in Pb-Pb collisions a further effect of the same relative magnitude 
which needs consideration is the distortion of the particle phase space by the Coulomb potential. This 
effect influences particles and antiparticles in opposite way, and has by factor two different strength for 
w-quark (charge +2/3|e|) and (d, s)-quarks (charge — l/3|e|). Because Coulomb-effect acts in opposite 
way on u and d quarks, its net impact on \ q is relatively small as we shall see. 

However, the Coulomb effect distorts significantly the expectation regarding \ s — > 1 for strangeness- 
deconfined source with vanishing net strangeness. The difference between strange and anti-strange 
quark numbers (net strangeness) allowing for a Coulomb potential within a relativistic Thomas-Fermi 
phase space occupancy model [j24|, allowing for finite temperature in QGP is: 



(N s - N s ) 



d?rd?p 



(3) 



which clearly cannot vanish for V in the limit A s — > 1. In Eq. (J3j) the subscript -Rf on the spatial 
integral reminds us that only the classically allowed region within the fireball is covered in the integration 
over the level density; E = y/m 2 + p 2 , and for a uniform charge distribution within a radius Rf of charge 

, 2" 



V 



3 Z f e 2 
2 R f 



Z f e 2 



1 



3 (flf) 



for r < R 



(4) 



for r > Rf 



One obtains a rather precise result for the range of parameters of interest to us (see below) using 
the Boltzmann approximation: 



(N s - N 3 ) = 7s 




-E/T 




\ — \-l - — 

A s e3T - A„ e st 



(5) 



The Boltzmann limit allows also to verify the signs: the Coulomb potential is negative for the negatively 
charged s-quarks with the magnitude of the charge, 1/3, made explicit in the potential terms in all 
expressions above. It turns out that there is always only one solution, with resulting A s > 1. The 
magnitude of the effect is quite significant: choosing Rf = 8fm, T = 140 MeV, m s = 200 MeV (value 
of 0.5 < 7 S < 2 is irrelevant) solution of Eq. (jH) for Zf = 150 yields A s = 1.10 (precisely: 1.0983, 1.10 
corresponds to Rf = 7.87 fm). This result is consistent with one of the scenarios we reported earlier for 
Pb-Pb collisions ||. Thus we are reassured that the experimental data is very likely consistent with 
deconfined quark source, and hence a detailed verification of this hypothesis is needed. The remarkable 
result we find is that experimental data is only consistent with this value X s = 1.10 ± 0.02, see table [l]. 
Thus as before for the lighter system S-Au/W/Pb 0, 15, 0] we are finding that the source of strange 



hadrons (up to Coulomb-asymmetry) is governed by a symmetric, and thus presumably deconfined 
strange quark phase space. 

Let us briefly explain how we obtain the statistical parameters shown in table [1} with E; t = 
l i + P 2 = Y m i + v\ cosh y we integrate over the transverse momentum range as given by the ex- 
periment, see table ||. To obtain the relative strengths of centrally produced particles we consider only 
central rapidity region y ~ 0. We allow all hadronic resonances to disintegrate in order to obtain the 
final relative multiplicity of 'stable' particles required to form the observed particle ratios. We show 
results of four main models denoted A, B, C, D in top section of table [l], arising describing the data 
shown in the first 4 columns of table ^. In this group we successively relax the chemical variables from 
their tacit values (= 1). In each step the number of degrees of freedom decreases by one, yet as described 
by x 2 / n do{ the confidence level becomes progressively better, and is indeed of impressive quality, with 
X 2 / n do{ = 0.3 when all 4 chemical variables are allowed to vary. We thus conclude that it is necessary 
in description of the particle abundance data to allow non-equilibrium abundances of light and strange 
quarks. 
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Table 1: Statistical parameters obtained seeking minimum of weighted least square difference with 
experimental results shown in table |2[ Values of A s in approach D s is the result of strangeness conser- 
vation constraint: asterisk * means a fixed (input) value, not a parameter. In case the temperature 
is fixed at the value obtained for S-Au/W/Pb collisions, and in case D p , the pressure in the hadronic 
phase space is fixed at 82 ± 6 MeV/fm 3 , the value obtained in S-Au/W/Pb collisions. In case F the 
four ratios with O-data are studied. 



Fit 


T f [MeV] 


A q 


A s 


7s 


7q 


Wdof 


X 2 /«dof 


A 


147 ± 3 


1.69 ± 0.03 


1* 


1* 


1* 


8 


16 


B 


142 ± 3 


1.70 ± 0.03 


1.10 ± 0.02 


1* 


1* 


7 


13 


C 


144 ± 4 


1.62 ± 0.03 


1.10 ± 0.02 


0.63 ± 0.04 


1* 


6 


4 


D 


134 ± 3 


1.62 ± 0.03 


1.10 ± 0.02 


1.27 ± 0.08 


1.84 ± 0.30 


5 


0.3 


D s 


133 ± 3 


1.63 ± 0.03 


1.09* ± 0.02 


1.98 ± 0.12 


2.75 ± 0.35 


5 


0.45 


D 4 


143* 


1.61 ± 0.03 


1.10 ± 0.02 


0.74 ± 0.06 


1.15 ± 0.18 


5 


0.8 


D p 


137 ± 4 


1.62 ± 0.03 


1.10 ± 0.02 


0.88 ± 0.07 


1.33 ± 0.10 


5 


0.5 


F 


334 ± 18 


1.61 ± 0.03 


1.12 ± 0.02 


0.09 ± 0.01 


0.18 ± 0.02 


8 2.4 



This is confirmed by the results we show in Fig. |l|, where for a given If , for case D (solid line, no O or 
strangeness conservation) and case D s (dashed line, no O, with strangeness conservation), with all the 
other parameters obtained finding minimum of weighted least square theory-experiment difference at 
given temperature. The locations of the best x 2 / Ti dof are indicated by vertical lines. The top sections of 
the figure shows that both chemical non-equilibrium parameters 7 S , 7 q > 1 in a wide range of freeze-out 
temperatures, indeed the values are slightly higher for the smaller Tf that are normally more favored on 
intuitive grounds (lower freeze-out particle density). However, freeze-out at Tf > 145 MeV would allow 
values 7s, 7 q < 1. It is reassuring that the analysis with matter flow IffiJ] even more clearly favors the 
high values of 7 S , 7 q > 1 we found here. We note also the constancy of the parameter A s (unconstrained 
solid line result) which assures us that the Coulomb effect we described cannot be ignored. We also 
note the counter-intuitive result that the energy per baryon at freeze-out obtained from the properties 
of the hadronic phase space using the statistical parameters, is dropping as the freeze-out temperature 
increases in the region of interest here. 

3 Hadronic particle abundances and phase space properties 

Since we have the statistical parameters in hand, we can as alluded to above, not only obtain the particle 
abundances within the kinematic cuts, shown in table H but also study the physical properties of the 
hadronic system at freeze-out, shown in table ||. 

We first observe, comparing the experimental data shown in left hand portion of table ^ with 
the different theoretical results presented on the right hand side, that we are more successful in the 
description of particles above the horizontal line. Below, in the bottom section of table ^ we group results 
comprising O and O-particles, which seem not to follow the same systematics, a fact already reported 
with respect of their spectral temperature by the WA97 collaboration A possible hypothesis is that 
a good fraction of these particles are made in processes that are different in nature than those leading 
to the other particle abundances, and hence we excluded O and O-particles from all but one of the 
approaches presented, denoted F. Moreover, if some strange particles hadronize separately, one cannot 
demand that the remaining particles balance strangeness exactly and hence we did not in general enforce 
strangeness conservation, except in case D s . That case is only slightly worse than D, which implies that 
the source for both Q and O is nearly symmetric with respect to abundance of strange and anti-strange 
quarks. This is consistent with the observation that much of the significant asymmetry in the ratio Q/O 
arises from the Coulomb effect we described above. We note that model F shown in table |], including 
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the four O-particle data points yields x 2 / n dof = 2.4 for radof = 8, the mathematical confidence level is a 
few percent, it can safely be assumed that our approach is not adequately accounting for the production 
of O-particle. 

We address extensively in our study the different constraints, as indicated by the subscript in all 
tables: 

1. D s includes the requirement of strangeness conservation, i.e., the hadronic phase space has to 
contain for the given statistical parameters as many s- as s-quarks. This is most conveniently ac- 
complished by finding the value of \ s which balances strangeness in terms of the other parameters, 
and thus, though not fitted, the value shown in table || displays an error, derived from the errors 
determined determining the other statistical variables. We note that the phase space occupancies 
change drastically between cases, D and D s , see table |], however the value 7 s /7 q changes from 
0.69 for D to 0.72 for D s . In actual numerical procedure we took advantage of this stability in 
7 s /7q-ratio, using it as a parameter. More generally, we note that all acceptable models shown in 
table H yield "f s /j q = 0.68 ± 0.05, which is consistent with the result of model C for 7 S , where the 
tacit assumption ^ q = 1 is made. 

2. There is considerable reason to seek a comparison of the Pb-Pb system with the analysis of 
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Table 2: Particle ratios: Experimental results for Pb-Pb, references and kinematic domain in first four 
columns, and different model results corresponding to those shown in table |] in the remaining columns; 
Asterisk * means the corresponding experimental result is not fitted. 



Ratio 


Ref. 


cuts[GcV] 


Exp. data 


Fit A 


Fit B 


Fit C 


Fit D 


Fit D s 


Fit Dt 


Fit Dp 


Fit F 


S/A 




LI 




p ± > 0.7 


0.099 ± 0.008 


0.130 


0.138 


0.093 


0.095 


0.098 


0.095 


0.093 


0.107 


H/A 


1 


LI 


1 


p± > 0.7 


0.203 ± 0.024 


0.361 


0.322 


0.198 


0.206 


0.215 


0.201 


0.201 


0.216 


A/A 


1 


L] 


1 


p±_ > 0.7 


0.124 ± 0.013 


0.126 


0.100 


0.121 


0.120 


0.119 


0.123 


0.120 


0.121 






" 




p± > 0.7 


0.255 ± 0.025 


0.351 


0.232 


0.258 


0.260 


0.263 


0.262 


0.258 


0.246 


(S+E) 

f \ \ W 

(A+A) 




3 




p± > 1. 


0.13 ± 0.03 
11.9 ± 1.5 


0.169 
6.42 


0.169 
6.3 


0.114 
10.4 


0.118 
9.89 


0.122 
9.69 


0.116 
10.1 


0.115 
10.2 


0.120 
16.1 


K + /K~ 




3 






1.80 ± 0.10 


2.27 


1.96 


1.75 


1.76 


1.73 


1.73 


1.77 


1.62 


P/P 




3 






18.1 ±4. 


20.5 


22.0 


17.1 


17.3 


17.9 


16.6 


17.3 


16.7 


A/p 




3 






3. ± 1. 


2.97 


3.02 


2.91 


2.68 


3.45 


2.91 


2.96 


0.65 


K°/B 










0.183 ± 0.027 


0.298 


0.305 


0.224 


0.194 


0.167 


0.211 


0.214 


0.242 


h~/B 




3 






1.83 ± 0.2 


1.30 


1.47 


1.59 


1.80 


1.86 


1.55 


1.72 


1.27 


n/E 




LI 




p± > 0.7 


0.192 ± 0.024 


0.115* 


0.119* 


0.080* 


0.078* 


0.080* 


0.082* 


0.078* 


0.192 


fi/5 


1 


L] 


1 


p± > 0.7 


0.27 ± 0.06 


0.32* 


0.28* 


0.17* 


0.17* 


0.18* 


0.18* 


0.17* 


0.40 


n/n 




~ 




p ± > 0.7 


0.38 ± 0.10 


1.00* 


0.55* 


0.56* 


0.57* 


0.59* 


0.56* 


0.56* 


0.51 


(n+n) 

(H+H) 


1 


ii 


I 


p_L > 0.7 


0.20 ± 0.03 


0.17* 


0.15* 


0.10* 


0.10* 


0.10* 


0.10* 


0.10* 


0.23 



S-Au/W/Pb reactions which we reported earlier Thus we consider model in which the 
freeze-out temperature is fixed at the value we found in S-Au/W/Pb reactions, Tf = 143 MeV. In 
model Dp, the pressure of the hadron phase space is chosen at the value we found in S-Au/W/Pb 
reactions, P = 82MeV/fm 3 . As judged by x 2 / n dof all Dj-models are possible, and the resulting 
particle multiplicities presented in the columns of table |2| differ only in minute detail. The two 
models and D p which test consistency with the smaller S-Au/W/Pb reactions are well within 
the allowable error. This consistency implies the possibility that the matter formed in these two 
very different systems hadronize in a rather similar fashion, though collective surface flows are 
very different. 

To resolve if there is universal freeze-out we have to consider the physical properties of the fireball. 
While the statistical parameters shown in table [l] can vary strongly from model to model, we find that 
the implicitly determined physical properties of the hadron source are more stable. In table |3| we show 
for the 8 models along with their temperature the specific energy and entropy content, and specific 
anti-strangeness content, along with specific strangeness asymmetry, and finally pressure evaluated by 
using the statistical parameters to characterize the hadronic particle phase space. We note that it is 
improper in general to refer to these properties as those of a 'hadronic gas' formed in nuclear collisions, 
as the particles considered may be emitted in sequence, and thus there never is a stage corresponding 
to a hadron gas phase. However, in the event such a stage exists, we also evaluated (see last column in 
table H) the volume of the hadron gas source at chemical decoupling. In order to obtain this extensive 
property, we used the net baryon number in the fireball being (B — B) = 372 ±10, as stated in [pLOf] . Note 
that a spherical source corresponding to the best model D would have a source radius 9.6 fm, which in 
turn can be checked to be exactly in agreement with deconfined strangeness conservation as described 
by Eq. (|3]), given the statistical parameters, and m s ~ 200 MeV. Other interesting conclusions arising 
in view of these results and shown in table || are: the specific energy content E{/B is well within the 
expectations based on the collision energy content per nucleon (8.6 GeV) and hence this result confirms 
firmly the hypothesis that the energy stopping and baryon number stopping in the fireball are very 
similar. The specific strangeness content of the Pb-Pb collision fireball is, by about 20% smaller than 
S-Au/W/Pb result. 
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Table 3: Physical properties of hadronic final state phase space (specific energy, entropy, anti- 
strangeness, net strangeness, pressure and volume) derived from statistical parameters shown in table 
|l| which describes the Pb-Pb data in table ||. Asterisk * means fixed input. 



rrt 


if [MeVJ 


bji/B 


bf/B 


Sf/B 


(Sf - Sf)/B 


Jr'f [CjreV/fm J J 


V{ [tm \ 


A 


147 ± 3 


6 6D + 40 


37.0 ± 3 


92 + 05 


29 + 02 


068 + 005 


6429 + 500 


B 


142 ± 3 


7.13 ± 0.50 


40.9 ± 3 


1.02 ± 0.05 


0.21 ± 0.02 


0.053 ± 0.005 


8994 ± 600 


C 


144 ± 4 


7.75 ± 0.50 


41.7 ± 3 


0.70 ± 0.05 


0.14 ± 0.02 


0.053 ± 0.005 


10242 ± 800 


D 


134 ± 3 


8.33 ± 0.50 


46.8 ± 3 


0.61 ± 0.04 


0.08 ± 0.01 


0.185 ± 0.012 


3619 ± 250 


D s 


133 ± 3 


8.72 ± 0.50 


48.1 ± 3 


0.51 ± 0.04 


0* 


0.687 ± 0.030 


1134 ± 100 


D t 


143* 


7.63 ± 0.45 


41.4 ± 3 


0.68 ± 0.05 


0.12 ± 0.01 


0.072 ± 0.005 


7517 ± 500 


D p 


137 ± 4 


8.05 ± 0.50 


44.7 ± 3 


0.67 ± 0.05 


0.13 ± 0.01 


0.082* 


7090 ± 500 


F 


334 ± 18 


9.79 ± 0.50 


24.1 ± 2 


0.78 ± 0.05 


0.06 ± 0.01 


1.64 ± 0.006 


2303 ± 250 



4 Conclusions 

We have presented detailed analysis of hadron abundances observed in central Pb-Pb interactions at 158 
A GeV in terms of thermal equilibrium and chemical non-equilibrium phase space model of (strange) 
hadronic particles. We assumed formation of a thermal dense matter fireball of a priori unknown 
structure, which explodes and disintegrates into the final state hadrons. We have presented several 
excellent descriptions of all abundance data which at present comprise 5 or more independent degrees of 
freedom, yielding a family of models with acceptable confidence level. The physical statistical parameters 
obtained characterize a strange particle source which, when allowing for Coulomb deformation of the 
strange and anti-strange quarks, is exactly symmetric, as is natural only for a deconfined state. While 
the statistical parameters shown in table [l] can vary widely there is no way to distinguish with naked 
eye the different approaches D, D s , D^, D p inspecting the particle abundances shown in table [2|. It is 
important to take note that along with X q = 1.62 db 0.03, A s = 1.10 ± 0.02 there also is a stable value 
ls/lq = 0.68 d= 0.05 under the different strategies one may follow to analyze the experimental data. The 
chemical freeze-out temperature allowing for the systematic uncertainty seen in the acceptable group of 
models in table |] is Tf = 138 ± 7 MeV, and this implies that the freeze-out baryochemical potential is 

= 199 ± 3 MeV. The error here is small, since the best values Tf, \ q are anti-correlated. 

Given these statistical parameters we have also evaluated the physical properties of the hadronic 
particle phase space, such as energy, entropy and baryon number. The results shown in table || describe 
the properties of the final state. These correspond nearly exactly to the initial state conditions, con- 
firming the consistency of our approach and validating the reaction picture applied. This part of our 
analysis confirms that the reaction proceeds by the way of the formation of a dense fireball comprising 
highly excited hadronic matter. In consistency with this we obtain values of \ s which exactly match 
expectations for strangeness balance in QGP, allowing for the Coulomb effect within the particle source 
of the size R{ = 9.6 ± 2fm. We compare conditions of the particle source for the two systems Pb-Pb 
and S-Au/W/Pb and find that both can be seen as hadronizing in same physical conditions. 

We have compared the lighter system S-Au/W/Pb with the current study of Pb-Pb, selecting 
comparable freeze-out conditions (models and D p ). We find, see table ||| and |Q] that the different 
physical properties of the two hadronic source agree. This, along with the strange phase space symmetry 
and the Coulomb effect, we believe that the sole possible interpretation the formation of a deconfined 
phase in the initial stages of the collision, which subsequently evolves and flows apart till it reaches 
the universal hadronization point, with many similar physical properties, independent of the collision 



system. System dependent will certainly be the surface collective velocity v c [21], however, our analysis 
was organized such that this vector field did not enter here in a significant way. 

We have begun, using quark-gluon plasma equations of state which incorporate the perturbative 
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corrections and thermal masses, to study detailed scenarios of QGP formation and evolution that leads 
to the freeze-out properties we obtained here. The important preliminary finding is that it is possible to 
find QGP-fireballs that naturally lead to the results obtained studying the experimental hadron abun- 
dance data, and thus the QGP hypothesis presented here is also consistent with our current theoretical 
understanding of the QGP equations of state. 
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